Entanglement and the Quantum Brachistochrone Problem 
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Entanglement is closely related to some fundamental features of the dynamics of composite quan- 
tum systems: quantum entanglement enhances the "speed" of evolution of certain quantum states, 
as measured by the time required to reach an orthogonal state. The concept of "speed" of quantum 
evolution constitutes an important ingredient in any attempt to determine the fundamental limits 
that basic physical laws impose on how fast a physical system can process or transmit information. 
Here we explore the relationship between entanglement and the speed of quantum evolution in the 
context of the quantum brachistochrone problem. Given an initial and a final state of a compos- 
ite system we consider the amount of entanglement associated with the brachistochrone evolution 
between those states, showing that entanglement is an essential resource to achieve the alluded 
time-optimal quantum evolution. 
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The celebrated brachistochrone problem has played a 
distinguished role in the history of mechanics [l[ . A quan- 
tum dynamical optimization problem akin to the classical 
brachistochrone one has attracted the attention of sev- 
eral researchers in recent years 01,11 HI- This problem 
deals with the Hamiltonian generating the optimal quan- 
tum evolution \if)(t)) (in the sense of the one requiring the 
shortest time r) between two prescribed states \ipi) and 
\ipp)- It has been recently established that the entangle- 
ment features of quantum states are important in connec- 
tion with dynamical optimization problems @, 0, H, Q . 
Previous research done on this subject has been focused 
upon the study of the time needed to reach an orthogonal 
state by systems evolving under a given Hamiltonian [Io| . 
The details of these kind of studies depend strongly on 
the particular Hamiltonian considered, and each case re- 
quires a separate treatment (see, for instance, d,Q). The 
aim of the present contribution is to revisit the connec- 
tion between entanglement and optimal evolution from 
the different (but related) point of view provided by the 
quantum brachistochrone problem. An important advan- 
tage of the brachistochrone approach is that it allows for 
a more general and unified investigation of the alluded 
connection, going beyond the separate analysis of indi- 
vidual cases. 



of dimensions Na and Nb, respectively). We want to 
quantify the amount of entanglement involved when im- 
plementing the evolution of an initial state to a fi- 
nal state \i/>f) m the shortest possible time. Following 
3 , we are going to consider the optimal evolution under 
the constraint that the difference between the maximum 
and minimum eigenenergies of the Hamiltonian generat- 
ing the unitary transformation \ipj) — > \4>f) — 
be less or equal to a given constant energy 2uj. This 
constraint is imposed due to the following reason: if the 
differences between the eigenenergies of the Hamiltonian 
are arbitrarily large, then it is easy to implement a quan- 
tum evolution connecting the alluded states and taking a 
time t that can be made arbitrarily small. The optimal 
time evolution is given by [3| 
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This expression satisfies \ip(0)) = and \ip(r)) = \ipi 
where 



I. ROLE OF ENTANGLEMENT IN 
TIME-OPTIMAL QUANTUM EVOLUTION 

As already mentioned, our aim is to explore the con- 
nection between entanglement and the quantum brachis- 
tochrone evolution of a bipartite system constituted by 
two subsystems A and B (with associated Hilbert spaces 



To determine the parameter 0, the final state \i/)p) has 
to be written in the form 

|<M = cos + e *(*+*/2) sin ^|^), (3) 
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where \tpj) is a state orthogonal to the initial state 
contained in the two-dimensional subspace (of the full 
Hilbert space) spanned by the initial and final states. 
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Since both \tpj) and \tpp) are specified, the values of <f> 
and 9 can be regarded as known, the latter being the 
angle of separation of these two states. We are going to 
focus upon the case where the initial and final states are 
orthogonal, thus having 9 = it and, 
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From the point of view of the physics of information and 
computation this case is the most interesting because the 
evolution of a computer devise into a state orthogonal 
to the initial one can be identified with an elementary 
information processing step [ll[. 

In order to assess how much entanglement is involved 
in the brachistochrone evolution we are going to compute 
the time-average of the entanglement £\ip{t)] during the 
optimal time evolution from \tpj) — > \iPf), 



£(t)dt. 



(5) 



The idea of time averaged entanglement has been recently 
discussed by several researchers and was found to consti- 
tute a useful concept to study a variety of problems (see, 
for instance, 12j and references therein). To facilitate the 
computations we use the linear entropy as entanglement 
measure, 
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where pa is the marginal density matrix associated with 
one of the two subsystems (the one with the Hilbcrt 
space of lower dimension), pa = Trs[p(t)], and the 
global density matrix describing the bi-partite system is 
p(t) = \ip(t))(ip(t)\. Since the entanglement £{t) depends 
upon time only through the quantity £ = x ' ^ ne integral 
giving the time-averaged entanglement can be put under 
the guise 
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implying that {£) is independent of oj and since r — ^ , 
this means that (£) is also independent of the absolute 
time taken. 

Highly asymmetric states can evolve optimally without 
entanglement. This is clearly illustrated by the extreme 
case corresponding to factorizable initial and final states 
of the form 



where one of the subsystems (in this case, the second 
subsystem) is in the same state |0o) at the beginning 
and at the end of the process. It is plain that in such 
circumstances the brachistochrone evolution can be im- 
plemented without entanglement. Indeed, the optimal 
evolution is given by the time dependent, separable state 



$(i)> = \<f>(t))®\<h), 
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\4>(t)) being the optimal evolution of the first subsystem 
connecting the states \4>i) and \<fe). It is also clear that 
in this highly asymmetric setting only one of the sub- 
systems (the first one) is evolving and, consequently, we 
are essentially dealing with the evolution of a single sys- 
tem. That is, the composite nature of the total system 
plays no role in an evolution like ([9]) . This is fully consis- 
tent with similar results reported by Giovannetti, Lloyd 
and Maccone in Q. There it was pointed out that, for 
composite a system with non-interacting subsystems, as- 
symetric non-entangled states in which all the energetic 
resources are concentrated in one subsystem can saturate 
the fundamental lower bound for the time required to 
reach a state orthogonal to the original one. On the con- 
trary, symmetric non-entangled states with evenly shared 
energetic resources do not saturate the alluded bound. 

States of the form f3|) are the only pairs of orthogo- 
nal initial and final states such that the time averaged 
entanglement of the concomitant optimal evolution van- 
ishes. This can be seen as follows. Suppose that the time 
averaged entanglement is zero. That implies that the en- 
tanglement is zero at all times t 6 [0, r]. In particular, 
the initial and the final states must be separable. Since 
they are also orthogonal we can assume, without loss of 
generality, that the initial and final states are, respec- 
tively, of the form \ipi) = |0) <g> \<j> r ) and \iPf) = |1) ® | </>«)■ 
Now, the time dependent state must be separable also at 
all intermediate times t € [0, r]. In particular, it must be 
separable at the time corresponding to ^ = \. That is, 
the state (1/V2)(|0) ® \<f> r ) + |1) (g) \(f) s )) has to be sep- 
arable. But for this state we have £ = 1 — \((f> r \(f>s)\ 2 ■ 
Consequently, the entanglement of this state will be zero 
only if |(0 r |(/> s )| = 1, and so the pair of initial and final 
states is of the form ([5]) . 



II. TWO-QUBIT SYSTEMS 

The features of the optimal evolution which are asso- 
ciated with the composite nature of the system (such as 
the role of entanglement) are more clearly seen when we 
study the evolution of symmetric states. For this reason 
we are now going to pay special attention on the optimal 
evolution of two-qubit systems associated with symmet- 
ric initial and final states. 

First we are are going to consider three specific cases 
of optimal evolution of symmetric sates of two qubit sys- 
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terns associated, respectively, with the following three 
pairs of initial and final states, 

• (i) |00) -> £^{|01) + |10)} + sina|ll) 

• (n) T5{|01) + |10)}-^{|00) + |11)} 

• R ^{|00)-i|ll)}^^{i|00)-|ll)}. 

Let us discuss each example separately: 
(i). In this case we have the evolution of a separable 
state into one of a family of possible final states parame- 
terized by the real parameter a. According to the value 
of this parameter we may have for the final state either 
another separable state (a — an intermediately en- 
tangled state (a € (0, ?)) or a maximally entangled state 
(a = 0). Let C = X' Then the time dependent, opti- 
mally evolving state, the concomitant entanglement, and 
its time average are respectively given by, 

cos cx 

|*(£,a)) = cos£|00)+sin£— [|01) + |10)] 
+ sin£ sin a |11), 
£(£, a) — ^cos(a) 2 sin(£) 2 — sin(a) sin(2£)^ , 

(£(a)) = ^ cos(q) 4 — — cos(a) 2 sin(a) 

8 7T 

+ ^shr(a) 2 . (10) 
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FIG. 1: Plot of (£(a)) as a function ofa,a£ [0, §]. 

Figure 1 is a plot of the time averaged entanglement 
(£(a)) as a function of the parameter a. For |00) — > |11) 
(a = ^) we have the optimal time evolution of two 
qubits from an initial symmetric separable state into 
a final symmetric separable state orthogonal to the 
initial one. In this case we have (£) = i, which is the 
maximum time-average entanglement in this family of 
evolutions. On the other hand |00) -> ^[|01) + |10)] 
(a = 0) constitutes an instance of the optimal evolu- 
tion of a separable state into a symmetric maximally 
entangled state. In this case we have that (£) = |. Of 



particular interest is the minimum of the time-average 
entanglement within the family, which happens when 

|00) -> -^[|01) + |10) + |11)] (a = arcsin^)) and 

gives {£) = 0.088298. 

(ii) . This constitutes a particular instance of the opti- 
mal transformation of a symmetric maximally entangled 
state to another symmetric maximally entangled state. 
The time dependent entanglement and its average are, 
respectively, 

£(i)=cos 2 (^), (£) = \. (11) 

In this particular case the same amount of time-averaged 
entanglement is needed to transform a maximally 
entangled state into another maximally entangled state 
as is needed to transform a separable state into another 
separable state. 

(iii) In this case the optimally evolving state is maximally 
entangled at all times, consequently giving £ = 1. 

III. TYPICAL ENTANGLEMENT PROPERTIES 
OF TIME-OPTIMAL EVOLUTIONS OF 
TWO-QUBIT SYSTEMS 

Even when considering only two qubits and restricting 
\ipi) and IV'f) to be symmetric and orthogonal, the ex- 
pression for (£} becomes quite involved, thus making an 
analytic treatment of certain properties of the optimal 
evolution unpractical. 

In order to explore the typical entanglement features of 
optimal evolutions in two-qubits systems we are going to 
generate random pairs of symmetric (orthogonal) initial 
and final states, 

= ci|00) + c 2 -^(|01) + |10))+c 3 |ll) I 
\^ F ) = di|00) + d 2 -^(|01) + |10))+d 3 |ll),(12) 

with Y,i |c 2 | - E« = 1 and £i adf = 0. We treat 
the case of symmetrical states separately because, as we 
shall presently see, in this case the necessity of entangle- 
ment to implement time-optimal dynamics is somewhat 
stronger than in the general case. However, an analysis 
similar to the one performed here for symmetrical states 
can be implemented for other particular families of states. 
We are going to consider the statistical distribution, and 
the mean and minimum values of the time averaged en- 
tanglement (£) associated with the time-optimal evolu- 
tions connecting the pairs of states and \i/jf)- Inter- 
preting the triplets (01,02,03) and (efi, cfejCfe) as vectors 
in a three dimensional Hilbert space, we are going to gen- 
erate the random states (|12[) by applying random 3x3 
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unitary matrices U to the vectors (1,0,0) and (0,1,0). 
The random matrices U are gene rated uniformly accord- 
ing to the Haar measure (see [13| and references therein) . 
The corresponding probability distribution for the time 
averaged entanglement (£) is depicted in Figure 2 (solid 
line). This distribution, describing the (fractional) vol- 
ume of state-space associated with different values of 
(£), has a lower cut-off at (£) m in = 0.03415330, mean- 
ing that, in order to implement a time-optimal evolution 
between symmetric two-qubit states a finite, minimum 
amount of entanglement is needed. Such an evolution 
cannot be implemented without entanglement. On the 
other hand, there are time-optimal evolutions between 
symmetric states exhibiting a maximum time-averaged 
entanglement (£} m ax — 1.0. The mean value (over all 
possible optimal evolutions connecting symmetric states 
of two qubits is (£) mea n = 0.5. 



Symmetric states 

General states — 
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FIG. 2: Plot of the probability density function of the time av- 
eraged entanglement {£) associated with symmetric and gen- 
eral two-qubit states. 



We have also performed a similar calculation for 
time-optimal evolutions connecting two general (that is, 
not necessarily symmetrical) normalized and orthogonal 
states of two qubits, = ai|00) + a 2 |01) + a 3 |10) + 
a 4 |ll> and \ip F ) = 6i|00> + 6 2 |01) + 6 3 |10) + 6 4 |11). The 
vectors (oi, a-x, 03, (Z4) and (61, 62, &3, ^4) were obtained by 
acting with random unitary matrices U (generated uni- 
formity according to the Haar measure) upon the vectors 
(1,0,0,0) and (0,1,0,0), respectively. The correspond- 
ing probability distribution for the (£)-values is shown 
in Figure 2 (dashed line). In this case there are optimal 
evolutions with arbitrarily small values of the time av- 
eraged entanglement. However, the probability density 
P((£)) approaches zero as (£) — > 0. Indeed, as we have 
already shown, the only optimal evolutions with (£) = 
are those of the form ijHJ) where one of the subsystems 
does not evolve. It can be clearly appreciated in Figure 
2 that, as a general trend, more entanglement is involved 
in optimal evolutions between symmetric states than the 
one involved in general optimal evolutions (see also Table 



IV. SYSTEMS OF HIGHER DIMENSIONALITY 

We have also considered the connection between entan- 
glement and brachistochronc quantum evolution for two- 
qutrit and three-qubit systems (as entanglement measure 
for three-qubit states we used the average bi-partite en- 
tanglement associated with the three bi-partitions of the 
system into a qubit and a two-qubit subsystem). The 
corresponding probability densities for the time averaged 
entanglement associated with brachistochrone evolutions 
connecting random orthogonal states are depicted in Fig- 
ure 3. It transpires from this Figure that, both for two- 
qutrit and three-qubit systems, a considerable amount of 
entanglement is involved in typical brachistochrone evo- 
lutions. 




0.2 0.4 0.6 0.8 1 



<E> 



3.5 - 3 qubits 

3 - 




0.2 0.4 0.6 0.8 1 



<E> 



FIG. 3: Plot of the probability density function of the time 
averaged entanglement associated with time optimal quantum 
evolutions between symmetric and general orthogonal states 
of two-qutrit and three-qubit states. 



The minimum, maximum, mean and likeliest values of 
the time averaged entanglement associated with brachis- 
tochrone evolutions between orthogonal symmetric (S) 
and general (G) states (for two-qubit, three-qubit and 
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two-qutrit systems) are given in Table |TJ 





2 qubits 


3 qubits 


2 qutrits 




S 


G 


S 


G 


s 


G 




0.0342 


0.0000 


0.0304 


0.0000 


0.0256 


0.0000 


E ma x 


1.0000 


1.0000 


1.0000 


1.0000 


0.9914 


0.9915 


Pj mea n 


0.5000 


0.4000 


0.6667 


0.6667 


0.6429 


0.5999 


^likeliest 


0.50 


0.41 


0.71 


0.71 


0.69 


0.64 



TABLE I: Minimum, maximum, mean and likeliest time- 
averaged entanglement associated with brachistochrone evo- 
lutions of systems of 2 qubits, 3 qubits, and 2 qutrits. 



V. CONCLUSIONS 

We have investigated the role of entanglement in time- 
optimal (brachistochrone) evolution of composite quan- 
tum systems. We proved that, except for trivial cases 
in which only one of the subsystems is actually evolv- 
ing, brachistochrone quantum evolution between orthogo- 
nal states cannot be implemented without entanglement. 
We studied in detail, for two-qubits, three-qubits and 
two-qutrits systems, the amount of entanglement (as 
measured by the time averaged entanglement) involved 
in typical brachistochrone evolutions connecting orthog- 
onal initial and final states. In all cases we found that a 
considerable amount of entanglement is needed in order 
to implement brachistochrone evolutions. The present 
approach to the study of the connection between entan- 
glement and quantum time-optimal evolution may be re- 
garded as "global" in the sense that (for a given system) 
it is not based upon the separate analysis of the the dy- 
namics generated by different possible Hamiltonians. If 
we consider, for instance, a two-qubits system, the global 
approach shows that, in general, entanglement is a neces- 
sary resource to implement optimal quantum evolutions. 
Previous approaches, on the contrary, would require a 
separate treatment of each different Hamiltonian of the 
system at hand (say, a two-qubit system) . 

The two alluded strategies, the "global" and the 
"Hamiltonian-specific" one, actually complement each 



other. The global approach, based upon the brachis- 
tochrone evolution, establishes in a general and unified 
way that there is a connection between entanglement 
and optimal quantum evolution between orthogonal pure 
states: considering at the same time (for a given system) 
all possible optimum quantum evolutions it is seen that 
most of them involve a considerable amount of entangle- 
ment. On the other hand, the Hamiltonian-specific treat- 
ment permits, for each specific quantum Hamiltonian, a 
more detailed analysis of the aforementioned connection 

Summing up, in the present work we have established 
a definite connection between entanglement and brachis- 
tochrone evolutions connecting pairs of pure orthogo- 
nal states. Arguably, time-optimal evolutions between 
pure and orthogonal states are (at least from the con- 
ceptual point of view) the most important ones, since 
the initial and the final state are fully distinguishable 
and the alluded evolution can be identified with one el- 
ementary information-processing step [Tlj |. However, it 
is possible to incorporate within the present approach 
brachistochrone evolutions connecting non-orthogonal 
pure states, since these evolutions are perfectly well de- 
fined @. Moreover, it would also be possible to explore 
the case of optimal time-evolution between mixed states, 
using the techniques advanced in [f|. We plan to inves- 
tigate these extensions of the ideas advanced here in a 
future contribution. 
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